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Abstract 

Two families of quasi exactly solvable 2x2 matrix Schrodinger operators 
are constructed. The first one is based on a polynomial matrix potential and 
depends on three parameters. The second is a one-parameter generalisation of 
the scalar Lame equation. The relationship between these operators and QES 
Hamiltonians already considered in the literature is pointed out. 

1 Introduction 

In a recent paper a classification of 2 x 2 matrix quasi exactly solvable (QES) 
Schrodinger operators in one spatial dimension is attempted. This problem was first 
addressed in and further developped in P and . Here we consider a suitable class 
of finite dimensional vector spaces of polynomials in a real variable and we construct 
two families of operators preserving sub-classes of these vector spaces. The first family 
is related to one of the cases treated in the second, which is not considered in 
generalizes an equation considered in The two corresponding QES equations 

respectively constitute "coupled channel" generalisations of the anharmonic and Lame 
QES scalar equations. 



Following the basic idea of QES operators @], [0 we consider the finite dimensional 
vector space of couples of polynomials of given degree n and m in a real variable x. 
We slightly generahze this vector space by setting 

/ 



V = P 



(1) 



V{n) 

where V{n) denotes the set of real polynomials of degree at most n in x while P is 
a fixed invertible 2x2 matrix operator; P can be interpreted as a change of basis 
in the vector space V{n) © V{m). With such an interpretation it is reasonnable to 
choose P of the form 

P= j (resp. P= 1^^^ ° j) (2) 

for n < m (resp. m < n). In this paper, we limit ourself to scalar operators Pu (or 
P21) of the form 

^ d d 

P12 = /tOTT + ^1 + l^2XTr + I^SX (3) 
OX OX 

where kj are constants. The vector space defined in Eq. (6) of [|l| can be set in the 
form (|I|) with Kq = 1, /ti,2,3 = 0. 



2 Polynomial potential 

We consider an operator of the form 

Hiy) = -£^^2 + M,iy) (4) 

where M^ly) is a 2 x 2 hermitian matrix whose entries are even polynomials of de- 
gree at most six in y. This operator is a natural generalisation of the famous QES 
anharmonic oscillator |^ to 2 x 2 matrix operator. After the standard "gauge 
transformation" of H{y) with a factor 

<P{y) = y' exp -{ yl/^ + Piy^} (5) 
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and the change of variable x = y'^, the equivalent operator H{x) can be computed : 

H{x) = r\x)H{y)<p{x)\y=^ (6) 

Then we pose the problem : what is the most general choice of Mq such that H{x) 
preserves a vector space of the form (|I]),(^,(|^)?. 

The following solutions are obtained after straighforward calculations (we exclude 
the case where M^ly) is diagonal since it corresponds to a direct sum of two scalar 
QES operators) : 

e = or 1 , n = m — 2 , ki = K2 = = . (7) 

The corresponding potentials Mq has the form 

M^iy) = {Apy + 8piP2y' + (4p? - 8mp2 + 2(1 - 2e)p2)y^}^2 

+ (8p2?/^ + 4]9i)o-3 - 8mp2Koai (8) 

where cri,(T3 are the Pauli matrices, P2,Pi,i^o are free real parameters and m is an 
integer. In particular, the non diagonal term is parametrized by an arbitrary constant 
which cannot be suppressed because of the y-dependent term proportional to a^. If 
the parameter e is choosen as an arbitrary real number, then the potential Mg has a 
supplementary term of the form e(e — 

When the parameters of the case 1 of Ref. are choosen so that the potential 
is a polynomial matrix (i.e. 02 = ao = 0, ai = 1, /?o = 1/2 or 3/2 in Eq. (34) 
of |]l|) the potential reduces to the matrix M^ly) above. The way of obtaining this 
result here is slightly different because the method starts from the natural vector 
space V{m) © V{n). The more elaborated QES operators obtained in can also be 
produced by our technique but this is not aimed in this note. 

Let us also point out that the "gauge factor U" considered in is limited to be 
a function of the variable x and, therefore, is not supposed to contain any derivative 
operator like our operator P (see (^). This explains that the QES polynomial po- 
tential (P) was not found in 0]. With the restriction that U is a function of x only, 
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these authors correctly reach the conclusion that Hamiltonian preserving a space like 
V{n) © V{m) with \n — m\> 1 are essentially diagonal, in contrast with the present 
operator related to the case n — m = 2. 

3 Application to N-body hamiltonians 

By using the idea of the previous section, a matrix version of the QES many-body 
problem of Ref. P] can be constructed. Let us consider the Calogero Hamiltonian 
(we note it Hcai) supplemented by a matrix- valued potential V* : 

H-H,.+v = lj:i-^+x^,]+Y.i^^+y' (9) 

^ j=i ^-^j j<:i -^i) 

Along with P we assume V* to depend only on the variable r 

N N 

r^Y.i^j-Y){x.-Y) , F^E^. > (10) 
and we look for eigenfunctions of the hamiltonian (^) of the form 

^(x)=^o(x) T'exp-i^T^+p^r''} <f){r) (11) 
where ipo denotes the ground state of the standard Calogero system : 

N 

M^) = [E\x^-x,\)Yexpi-Xy2) , X'^J2^' (12) 

i<j j=l 

while (/)(r) represents a couple of polynomials in r. 

After a standard algebra, the operator acting on 0(r) can be isolated : 

h^r-^^ + {4r + 2b)^ + V* (13) 

and it can be shown that this operator preserves the space ([^) (again with n = m — 2, 
Ki = K,2 = 1^3 = 0) provided V* is of the form 

V*{'r) = -Plr^ + 2p2(l - + {a- 2p2(X3)r + (1 - pO^s + ^ + 2mKo(Ti (14) 

r 
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with the definitions 

a = pi(2-pi)+p2(2m + 3e-l + 6) , b = ^il + uN){N -1) , 7 = 2e(e-l + 6) . (15) 

As a consequence, (|^),(0) constitutes a QES matrix extension (labelled by the pa- 
rameters pi,P2, e) of the exactly solvable Calogero hamiltonian. 



4 Lame type potential. 



As a second example, we consider the family of operators 
rf2 



Hiz) 



+ 



Ak'^sn^ + 5(1 + A;2) /2 2^A;cn dn 

2^A;cn dn Ck\n^ - 5(1 + k^) /2 



(16) 



where A, C, 5, ^ are constants while sn, cn, dn respectively abbreviate the Jacobi ellip- 
tic functions of argument z and modulus k |l^ 



sn(2;, k) , cn(2;, k) , dn(z, /c) . 



(17) 



These functions are periodic with period 4K{k),4K{k),2K{k) respectively {K{k) is 
the complete elliptic integral of the first type). The above hamiltonian is therefore 
to be considered on the Hilbert space of periodic functions on [0,AK{k)]. For com- 
pleteness, we mention the properties of the Jacobi functions which are needed in the 
calculations 



cn^ + sn^ = 1 



dn^ + A;^sn^ = 1 



-sn = cn dn , — cn = — sn dn , — dn = —k sn cn 



(18) 
(19) 



dz dz dz 

The relevant change of variable which eliminates the transcendental functions sn. 



cn, dn from (|l^) in favor of algebraic expressions is (for k fixed) 

X = sn^{z, k) 



(20) 



In particular the second derivative term in (|16[) becomes 



dz^' 



rf2 

4x(l - x)(l - k'^x)— + 2(3A;2x2 - 2(1 + k'^)x + 1 



rf2 



X 



dx 



(21) 



Several possibilities of extracting prefactors then lead to equivalent forms of (jlED, 
say H{x), which are matrix operators build with the derivative d/dx and polynomial 
coefficients in x. The requirement that H{x) preserves a space of the form (|lD leads 
to two possible sets of values for A, C, 9 (we do not consider the case 9 = since it 
corresponds to two decoupled scalar Lame equations). 

Case 1 

A = Am^ + 6m + 3 - 5 
C = Am^ + 6m + 3 + 5 
^ = i[(4m + 3)2-52]i 

The parameter 6 remains free, and also k which fixes the period of the potential. Four 
invariant spaces are available. In order to present them we conveniently define 

4m + 3 - (5 



Ri 



Am + 3 + 6'' 



(22) 



We have then 
Vi = 

V2 = 

Vs = 

V4 = 




k^Ri 



k^Ri 



Ri/k' 



(23) 
(24) 
(25) 
(26) 



Case 2 

A = Am^ + 2m + 1 - 5 

C = Arri^ + 2m + 1 + 5 
e = i[(4m+l)2-52]i 

The associated invariant vector spaces read, defining i?2 = (4m +1 — 5)/ {Am + 1 + 5), 
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V7 




(27) 
(28) 
(29) 



kyR2 (30) 



The operator (|T6|) was studied in for 5 = 1. For this particular value of S, the 

corresponding eigenvalue equation Hip = u'^ip determines the normal modes of the 
sphaleron classical solution in the Abelian Higgs model in 1+1 dimension. It 
therefore plays a crucial role in the understanding of the instabilities of the sphaleron 
in this model. 

The above results demonstrate that the remarkable algebraic properties of the 
Lame equation |jTO[ also hold for the operator (0), irrespectively of the value of 5. 
The associated eigenvalue equation therefore constitutes a (one parameter) 2x2 
matrix equation analog of the scalar Lame equation. 



5 Generalization 

The kind of operators presented in Sect. 3 can be generalized to matrix potentials of 
the form 



Hiz) 



e sn"i cn"2 dn°3 



6 sn"i cn"2 dn"^ l^fsn^l 



(31) 



where Vi, V2 are polynomials, ^ is a constant and aj are non- negative integers. 
The similarity transformation 



H[x) = U-\z)H{z)U{z) , U{z) = diag(sn^i cn^^ dn^^ sn^i cn^^ dn^^) (32) 
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sets the operator (^Tj) into a form with polynomial coefficients in the variable x = sn^ 
provided 

• Pj^lj = or 1 , j = 1,2,3 

• aj ± {j3j — 7j) = non-negative even integer , j = 1, 2, 3. 

After making a choice of aj,Pj,jj satisfying the above conditions, the possible 
forms of Vi, V2 and of P, m, n in Eq. have to be determined in order for H[z) to 
be QES. 

Taking k = 0, the standard trigonometric functions are recovered : 



sn(z,0) = sinz , cn(z,0) = cosz , dn(z, 0) = 1 



(33) 



The periodic potential below, which is exactly solvable ||12|, furnishes a particular 
example of this type 



V{z) 



cos^(z) 



cos[z) sm[z) 

,2/ 



(34) 



cos{z) sm{z) sin (z) 

It determines the normal modes about some static solutions of the Goldstone model 
in 1+1 dimensions 



6 Concluding remarks 

The examples of operators presented above give evidences of the difficulty to classify 
the coupled-channel (or matrix) QES Schrodinger equations. The way of constructing 
the QES potential Mq in Sect. 2 further provides a clear link between the approaches 
|]l]] and to this mathematical problem; we hope that this note will motivate further 
investigations of it. 
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